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. Recently H.Georgi has introduced the concept of unparticles in order to describe the low energy 

physics of a nontrivial scale invariant sector of an effective theory. We investigate its physical effects 
;_h on the neutrino flux to be detected in a kilometer cubic neutrino telescope such as IceCube. We study 

the effects, on different observables, of the survival neutrino flux after through the Earth and the 
regeneration originated in the neutral currents. We calculate the contribution of unparticle physics 
to the neutrino-nucleon interaction and, then, to the observables in order to evaluate detectable 
, effects in IceCUbe. Our results are compared with the bounds obtained by other non-underground 

experiments. Finally, the results are presented as an exclusion plot in the relevant parameters of 
the new physics stuff. 
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_ ■ I. INTRODUCTION 
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f-*- ' The Standard Model (SM) for the elementary particles interactions has been successfully tested at the level of 
quantum corrections. In particular high precision and collider experiments have tested the model and have placed 
the border line for physics effects at energies of the order of ITeV [![. On the other hand, new physics effect in the 
neutrino sector have recently received an important amount of experimental information coming from flavor oscillation 
This fact is the first evidence of neutrino masses different from zero, and hence, of physics beyond the SM. In this 
way, the neutrino sector and in particular neutrino-nucleon interactions, could be the place where new physics may 
become manifest again. Although the SM has been successful to describe the world at short distances, as a low energy 
effective theory of phenomena at higher scales, it leaves several open questions, e.g.: it does not predict the number 
of families and the fermions masses, has several free parameters, the mass generation mechanism through the Higgs 
boson, where its mass is not predicted, is untested and still leaves open the hierarchy problem. In these conditions, it 
is believed that we should have some kind of physics beyond the SM, which is called New Physics (NP)[l[. The search 
of NP proceeds mainly through the comparison of data with the SM predictions. The experimental way to look for 
NP effects in a model independent fashion is to construct observables that can be affected by this new physics and 
then compare the measurements with the mentioned SM expectation. Certain types of NP can already be present at 
the TeV scale and could participate in neutrino-nucleon interactions. Hence, these NP effects could possibly become 
apparent in neutrino telescopes. 

In the other hand it is well-known that scale invariance has been a powerful tool in several branches of physics and 
the possibility of a scale invariant weak interacting sector with the low energy particle spectrum is not ruled out. In 
particular, H. Georgi 0] has proposed that a scale invariant sector which does not imply conformal invariance 0] with 
a non-trivial IR fixed point and coupled to the SM fields through the exchange of particles with a high mass scale Mjj 
may appear much above the TeV energy scale. Below this energy scale, this sector induces unparticle operators Ou 
with a non-integral scale dimension d u that in turn have a mass spectrum which looks like a d u number of massless 
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FIG. 1: Diagrams contributing to the neutrino-nucleon cross section. 



particles. The couplings of these unparticles to the SM fields and in particular to standard neutrinos (i.e. massless 
and left-handed particles) and quarks are described by the effective Lagrangian Q 

A k+l-du _ 

£e ff = m hviCv + C Al5 )f C%, (1) 

where A is the energy scale at which scale invariance emerges, the dimensionlcss coefficients Cy(A) are °f order 1 for 
neutrinos and quarks and / is the generic fcrmion spinor. The operators with lowest possible dimension have the most 
important effect in the low energy effective theory regime. In Eq.(l) we have only included the vector unparticles 
operators that couple with left neutrinos and both left and right quarks. Note that the left couplings to the 
neutrinos and quarks are taken equal. 

In the present work we follow the Georgi's approach where the Feynman propagators of the unparticle operator 0£ 
is determined by the scale invariance Q , 
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The scale dimension d u is restricted in the range 1 < d u < 2. Here, the condition d u > 1 is due to the non-intcgrablc 
singularities in the decay rate [j| while d u < 2 is due to the convergence of the integrals @ . 

II. THE CROSS SECTION NEUTRINO-NUCLEON AND THE UNPARTICLES CONTRIBUTION. 

In this section we consider the effective operator given in Eq.([T]) and calculate their contribution to the neutrino- 
nucleon inclusive cross section: 

vN — » fi + anything, (4) 

where N = — - — is an isoscalar nucleon. The corresponding process is pictured in Fig. [1] which has the SM charged 

and neutral current diagrams and the unparticle contribution. For charged currents the calculation is standard and 
we use it to compare our results with Q. For neutral currents we have included the contributions of unparticles. The 
corresponding coupling and propagator are given in equations (fT]) and J2]) respectively. 

The SM results for the scattering amplitude for charged currents muon-neutrino scattering is 



and the corresponding differential cross section reads 
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where for an isoscalar target we have the quark distribution functions 

n cc ( n 2s _ u v {x,Q 2 ) + d v (x,Q 2 ) u s (x,Q 2 ) + d s (x,Q 2 ) 

+ s s (x,Q 2 ) + b s (x,Q 2 ), 

QCC{ ^ Q2) = U s (x,Q 2 )+ds(x,Q 2 ) + Cs(x> Q2) + ts{x> Q 2 } 

Similarly, the SM neutral current amplitude is 

■ 2 

M ™ = ~ 9,2 (7WTm*\ V ^ PlV E ^(SIlPl + 9WQi, (8) 

zt WW ~*~ m Z) i=U.D 

where cw = cos9 Wl x-w = sin 2 Q w , gY = 1/2 — 2x\y/3, g£ = — 1/2 + x w /3, = —2xw/3 and gj^ = xw/3. 

From the effective interaction Eq.(JTJ) between unparticle and the SM fields we obtain the following four fermion 
amplitude 



where the left (cl) and right (cr) coupling constants are expressed in terms of the vector and axial vector coupling 
constants: 



cr = cy + ca 
cl — cy - ca, 



(10) 



where 



CV{A) = M^M^~ du ■ 
The total neutral current contribution, including the unparticle contribution, can be written as 
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where 
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The neutral current differential cross section is then 
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FIG. 2: Total cross section for the SM and for different values of the unparticles dimension d u and cl = cr = 0.01. 

where the corresponding parton distributions for a isoscalar target read 
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In Fig. [2] we show the behavior of the total cross section (cr t (£') = a cc (E) + a NC (E)) with the neutrino energy for 
different values of d u and cl — cr = 0.01. We can appreciate a considerable disagreement with the SM predictions, 
due to the unparticle propagator, particularly for low values of d u and low neutrino energy. This very disparate 
behavior do not directly translate to the neutrino flux due to strong regeneration effects, as we will see in the next 
section. 



III. THE SURVIVING NEUTRINO FLUX 

The surviving flux of neutrinos of energy E, with inclination 6 with respect to nadir direction, <&{E,ff), is the 
solution of the complete transport equation Q : 



d\jx$(E,T') 
dr' 



-o-*(J5) + dE 
Je 



, §(E',t') da NC 
E ™~ ${E,t') dE ' 



(16) 



where the first term correspond to absorption effects and the second one to the regeneration. Here, r = r(6) is the 
number of nucleons per unit area in the neutrino path through the Earth, 



/•2_Re cos 6 

t{6) = N A / p(z)dz, (17) 

Jo 

In order to find a solution for this equation we make the following approximation Q: we replace the fluxes ratio 
inside the integral of the second member by the ratio of fluxes that solve the homogeneous equation (i.e., only 
considering absorption effects) 
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FIG. 3: Upper: 5 in the SM. We also include the factors ^ a bs and ^ r 
d u = 1.3 and cl = cr = 0.01. 



Lower: 5* including the unparticle contribution for 



where 

A(E',E) = [a t (E')-a t (E)]T' 
Thus, integrating the transport equation we have 
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where 
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&o(E) is the initial neutrino flux considered as isotropic. Na is the Avogradro number, i?E is the radius of the 
Earth, 9 is the nadir angle taken from the down- goin g normal to the neutrino telescope and the earth density p(r) 
is given by the preliminary reference earth model [13| . It is important to mention that the solution of the transport 
equation, Eq. (|20[) is the first, but quite accurate, approximation of the iterative method showed in Ref . fl0| . 

In order to illustrate the general behavior of the solution we show in Fig. [3l for 9 = 0, the factor S 

S = = VaUE, 6)* reg (E, 9), (23) 

with 

* abs (E,9) = e -^ E ^ 

y reg (E,9) = e^(E,e)T(e) I ) 

In this figure we show S for the SM and for the unparticles contribution with d u = 1.3 and = cr — 0.01 and we 
have explicitly included both factors (\& a i« and & r eg) to see the compensation between absorption and regeneration 
effects. 

For d u > 1.3 the absorption and the regeneration practically compensate each other and, then, the values of S are 
very near to the SM values (Fig. [3] Upper). For d u < 1.3 the absorption and regeneration effects do not cancel that 
efficiently and then we have values of S that clearly differ from the SM value. In Fig. 2] we show S for different values 
of the unparticle parameters and for two angles between the directions of the neutrino beam and the nadir. 



IV. THE OBSERVABLE a(E) 

The angle a(E) and the related ratio rj(E) introduced in Ref.[ll[ are the observable that we shall use in this paper 
in order to study the impact of the unparticle physics on neutrino detection in a neutrino telescope such as IceCube. 
By definition a(E) is the angle that divides the Earth into two homo-event sectors. When neutrinos traverse the 
planet in their journey to the detector, they find different matter densities, and then, different number of nucleons to 
interact with. In this conditions, the number of neutrinos that finally arrive to the detector depends on the arrival 
directions, indicated by the angle 9 with respect to the nadir direction. If we consider only upward-going neutrinos 
of a given energy E, that is, the ones with arrival directions 9 such that < 9 < tt/2, there will always exist an angle 
a(E) such that the number of events for < 9 < a(E) equals that for a(E) < 9 < tt/2. 

Clearly, the value of a(E) is energy dependent. For low energies, the cross section decreases and the Earth becomes 
transparent to neutrinos. In this case a(E) — > tt/3 for a diffuse isotropic flux since this angle divides the hemisphere 
into two sectors with the same solid angle. Obviously for extremely high energies, where most neutrinos are absorbed, 
a{E) — > 7r/2, and for intermediate energies a(E) varies accordingly between these limiting behaviors. 

In order to define a(E) wc consider the expected number of events (muon tracks though charged currents v^N 
interactions) at IceCube in the energy interval AE and in the angular interval A9 that can be estimated as 

Af = n T T [ [ dndE„a cc (E)<$>(E,9), (25) 

J A6 J AE 

where tit is the number of target nucleons in the effective detection volume, T is the running time, and a cc (E) is 
the charged neutrino-nucleon cross section. We take the detection volume for the events equal to the instrumented 
volume for IceCube, which is roughly 1 km 3 and corresponds to tit — 6 x 10 38 . The function $>(E,8) in Eq. (|2"5|) is 
the survival flux which is the solution Eq. (|2"D)) of the complete transport equation Q. 

The definition of a(E) is essentially the equality between two number of events, thus, to a good approximation, 
for each energy bin all the previous factors cancel except the integrated fluxes at each side. In this way, a(E) can be 
defined by the equation 

ra(E) pir/2 

/ d8sm9e~< 7 '" {E)T{8) = / dd sin0 e - CT ^^ T W, (26) 

JO Ja(E) 

which is numerically solved to give the results shown in the Fig. [5] There we show the SM prediction for a(E) and 
the unparticles contribution for different values of the dimension d u . 

The main characteristics of a(E) have been reported recently in Ref.[ll|. It is worth to notice that a(E) is weakly 
dependent on the initial flux but, at the same time it is strongly dependent on the neutrino nuclcon cross-section. 
Hence, the use of the observable a(E) reduces the effects of the experimental systematics and initial flux dependence. 




FIG. 4: S for the SM and for different values of d u at two different angles with respect to nadir. 



Since the functional form of a(E) sharply depends on the interaction cross section neutrino-nucleon, if physics beyond 
the SM operates at these high energies, it will become manifest directly onto the function a(E). 

In order to evaluate the impact of the observable a(E) to bound new physics effects, we have estimated the 
corresponding uncertainties on the SM prediction for a(E). Considering the number of events as distributed according 
to a Poisson distribution the uncertainty can be propagated onto the angle cxsm(E). The number of events TV as a 
function of q;sm is 

TV = 2Tm T TAEa cc (E)<i>o(E) / d0sin0e -<Tc " (£;)T(8) , (27) 

Jo 

where we have considered the effective volume for contained events so that an accurate and simultaneous determination 
of the muon energy and shower energy is possible. For IceCube, it corresponds to the instrumented volume, roughly 
1 km 3 , implying a number of target nucleons nx — 6 x 10 38 . We have considered an integration time T = 15 yr which 
is the expected lifetime of the experiment. To propagate the error on TV to obtain the one on a, we note that 
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FIG. 5: The predictions for a(E) obtained for different values on d u 



and dividing by N we obtain 
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where for Poisson distributed events we have 



5N = VN 



(29) 
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In order to evaluate the errors on a(E), it is necessary to consider a level of initial flux $>q(E). Here we have added 
together the cosmological diffuse isotropic flux and the atmospheric one(see Fig. [7|). For the atmospheric flux, we 
have considered the one given in Ref.[lJ]. As for the cosmological diffuse flux, the usual benchmark is the so-called 
Waxman-Bahcall (WB) flux for each flavor, E^4>y^ B ~ 2.4 x 10~ 8 GcV cm _2 s _1 sr _1 , which is derived assuming that 
neutrinos come from transparent cosmic ray sources fl5| . and that there is an adequate transfer of energy to pions 
following pp collisions. However, one should keep in mind that if there are in fact hidden sources which are opaque to 
ultra-high energy cosmic rays, then the expected neutrino flux will be higher. 

On the other hand, we have the experimental bounds set by AMANDA. A summary of these bounds can be found 
in Refs. (T5. and as a representative value we take Ef, ~ 8 x 10~ 8 GcV cm~ 2 s _1 sr~ 1 . With the intention 

to estimate the number of events, we have considered an intermediate flux (INT) level slightly below the present 
experimental bound by AMANDA, 



E, 
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(31) 



As it was discussed in Ref.[Tl[ the interval for maximum sensitivity for a is 10 5 GeV < E < 10 7 GeV. However, as 
for lower energies the atmospheric flux grows and then the errors fall, we have considered as an energy window for the 
fits the interval: 10 3 GeV < E < 10 7 GcV. In Fig. [fj]we show our results for the observable a and the corresponding 
errors within the mentioned energy window. In Fig. [7] we show the used flux. 

In the same context, we can define another observable related to a(E). We consider the hemisphere < 8 < ir/2 
divided into two regions by the angle cxsm(E), TZi for < 9 < asAt(E) and 7?-2 for c<sm(E) < 6 < n/2. We then 
calculate the ratio r)(E) between the number of events for each region, 



No- 



(32) 



where Ni is the number of events in the region TZi and N2 the number of events in the region 72.2 • By using rj(E) 
the effects of experimental systematic and initial flux dependence are reduced. If there is only SM physics, then we 
have that the ratio tism(E) = 1. In order to estimate the capability of r]{E) to bound unparticle effects, we have 
considered the values of rj(E) along with their error bars in Fig. [8] as if they had been obtained from experimental 
measurements for i](E). We proceed, then, to perform a x 2 -a n alysis taking as free parameters the dimension d u and 
the constant c = cl = cr and considering as experimental point the SM values for rj(E) for the same energy bin used 
in Fig. [fJJ We define the x 2 function in the usual way, 
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FIG. 6: The prediction for asAi(E) and the statistical errors. 



where according to the definition of r](E) fEq. ([3"2")) ) the statistical errors are given by Sx(Ei) = 2/\fNl for events 
distributed according to a Poisson distribution. The function \ 2 is minimized to obtain the allowed region in the 
(d u , c) plane, which corresponds to the region below the curve A shown in Fig. [9j In the same figure we also include 
other bounds obtained from different processes. The curve B corresponds to bounds obtained from atomic parity 
violation (l8j. the curve C corresponds to bounds from the muon anomaly and curve E comes from low energy 
v e — e scattering [2(|. This last bound is significatively better than all the other bounds. We would like to stress 
that the curve E corresponds to v e — e scattering but we arc considering muon neutrinos, i.e., a different neutrino 
flavor. Thus, if we consider the unparticle interactions as non-universal then this constraint would not apply and the 
neutrino telescope bounds may be competitive. However, there are experimental data for the ratio between — N 
Neutral Current to Charged Current interaction R [2l[ : 



fi(gj= (34) 



In order to compare with the experimental value 



R = 0.320 ±0.010, (35) 

which is an average measurement for energies in the range 20 GeV - 200 GeV, we calculated the corresponding 
average value for R: 

{R} = mwL v RiE)dE m 

and compared it with the previous experimental value fEq. (|35p ) for different values of the unparticles parameters 
(d u , c). The allowed region corresponds to the points below the curve D of Fig. [5] This curve shows that the bound 
coming from R is stronger than the one IceCube could possible set in the future. 
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FIG. 7: The utilized flux obtained by adding the atmospheric and the isotropical cosmic flux. 

V. CONCLUSIONS 

In the present work we studied the effects of unparticles contributions to the neutrino-nucleon cross section on the 
survival neutrino flux in a neutrino telescope like IceCube. To do it, we considered an effective interaction between 
standard particles and unparticles. We have found a considerable disagreement with the SM prediction for the neutrino 
observables defined above, particulary for low values of d u and low neutrino energy. For moderate values of d u this 
disagreement tends to disappear due to a strong cancellation between absorption and regeneration. 

We have also studied the possibility to bound effects of unparticles contributions to the interactions between muon 
neutrinos and the nucleons of the Earth using the observable rj(E). In this context, we fitted the theoretical expression 
for rj(E) as a function of the d u and c taking as experimental data the SM values obtained for rj (risAi(E) = 1) along 
with the errors derived for a number of events distributed according to a Poisson distribution. The results are shown 
in Fig. [HI as a allowed region plot. Finally, in the same figure we have compared the limits that have been obtained 
for different authors for low energy processes. The bounds coming from the ratio R of neutral to charged currents 
(curve D) is more restrictive than the bound IceCube could set (curve A) and then new effects that IceCube may 
detect could not be attributable to unparticle physics. 
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